A layering model for superconductivity in the Borocarbides 
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We propose a superlattice model to describe superconductivity in layered materials, such as the 
borocarbide families with the chemical formula -Rr2B2C and RTHC, with R being (essentially) 
a rare earth, and T a transition metal. We assume a single band in which electrons feel a local 
attractive interaction (negative Hubbard-?7) on sites representing the TB layers, while f/ = on 
sites representing the RC layers; the multi-band structure is taken into account minimally through 
a band offset e. The one-dimensional model is studied numerically through the calculation of the 
charge gap, the Drude weight, and of the pairing correlation function. A comparison with the 
available information on the nature of the electronic ground state (metallic or superconducting) 
indicates that the model provides a systematic parametrization of the whole borocarbide family. 
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The study of Nickel borocarbides (BCs), with the gen- 
eral formula i?Ni2B2C, has eiven rise to a wealth of in- 
teresting information 0, S B ■ Superconductivity has 
been found for R = Sc, Y, Ce, Dy, Ho, Er, Tm, Lu, 
and Th, but not for R = La, Pr, Nd, Sm, Yb, U, Gd, 
and Tb. In addition to the 'high' critical temperature for 
some compounds - Tc > 15 K for i? = Sc, Y and Lu -, 
the other aspect which has drawn attention to these ma- 
terials is the coexistence, for R — Dy, Ho, Er, and Tm, of 
superconductivity and magnetic ordering of the 4/ local- 
ized moments. Structurally, the BCs are characterized by 
a stacking of alternating RC sheets and Ni2B2 layers Isj, 
with both its band structure 0, and the near isotropic 
resistivity being consistent with three-dimensional be- 
haviour. 

From the microscopic point of view, superconductiv- 
ity in these materials is mediated by phonons, as evi- 
denced by specific heat ^ and isotope effect ^ mea- 
surements. It would then appear natural to relate the 
suppression of superconductivity, as R and the transi- 
tion metal (T) are varied, in terms of the BCS param- 
eters, Qd, ]^{Ef), and V (respectively, the Debye tem- 
perature, the density of states at the Fermi level, and 
some measure of the electron-phonon coupling strength): 
ksTc — exp[— l/Af(i?F)V']. Since resistivity measure- 
ments indicate that V does not vary much with R in 
_RNi2B2C one is left primarily with Qd and N{Ef). 
While Qd generally increases as R goes from Lu to La, 
measurements of the Sommerfeld coefficient, 7 oc N(E-p), 
lead to 7lu ~ '^Iha. in the Ni series. However, as the 
transition metal is varied, this simple parametrization 
no longer accounts for the trend of experimental data in 
an unambiguous way. Indeed, LaPt2B2C is a supercon- 
ductor [and, for this Pt series, so are the compounds with 
R = Pr, Y l^, and, possibly [13, Nd], even though it has 
a smaller N{Ef) than non-superconducting LaNi2B2C: 
7 = 5 and 8 mJ/mol K^, respectively Another im- 
portant piece of information is provided by studies of a 
closely related phase, i?NiBC |llil[l3. The structural 
difference between i?Ni2B2C and i?NiBC is the number 
of i? C layers between the Ni2B2 stacks: one in the former 
and two in the latter. Since these 'double-layered' ma- 
terials have so far failed to display superconductivity in 



both magnetic and non-magnetic systems alike, the lay- 
ering structure seems to play a crucial role in the BCs. 

Considerable insight should therefore be gained by de- 
scribing the i?r2B2C phenomenology in terms of zero- 
temperature transitions between superconducting and 
normal ground states as R and T are varied. With this 
in mind, here we propose a model based on the expec- 
tation that superconductivity can only survive the huge 
4/ local moment (~ lO/ie) in, say, DyNi2B2C if pairing 
of conduction electrons occurs solely in the NiB layers. 
Accordingly, we adopt a simplified description of single- 
band electrons feeling an attractive interaction only when 
occupying the same site on a T2B2 layer; electrons on 
sites of the RC layers do not interact with each other. 
Our purpose here is to single out effects such as chemical 
composition and layering in the suppression of supercon- 
ductivity; therefore, at this stage we will not take into 
account effects arising from the 4/ electrons of R, such 
as possible localized magnetic moments and hybridiza- 
tion. The tetragonal structure suggests a further sim- 
plification, namely, to consider a one-dimensional super- 
lattice, an approach which has been successfully used in 
interpreting some magnetic properties of the BCs; in do- 
ing this, we still capture the essential features of the in- 
terplay between layering alon g th e c-axis, and attractive 
electron-electron interactions And, finally, we incor- 
porate minimally the presence of actually different bands 
by including a band offset between the layers; that is, we 
allow for the occurrence of some charge transfer between 
the layers. As we will see, the model provides an effec- 
tive parametrization, through which one is able to sys- 
tematize the data and predict which materials are likely 
to be superconducting. 

The above picture is realized by a one-dimensional su- 
perlattice generalization of the attractive Hubbard model 

Ellii, 



n = Y^ E,n^a - t (cLc^+i^ + B.C. j + ^ C/, 



(1) 



where the notation is standard, apart from considering 
position-dependent site and correlation energy profiles: 
Ui = —\U\ and Si — e for all La sites of the attractive 
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Table I: Criteria to determine the nature of the ground state 
in the attractive Hubbard superlattice, from the behavior of 
the charge gap, A, and of the Drude weight, D. 



Nature of the state 


A [Eq. ©] 


D [Eq. ©] 


Superconducting 






Metalhc 







Insulating 








(i.e., Ni2B2) 'layer', and Ui = et = Q for all Lq sites on 
the free (i?C) layer; the basis with N}, = La + Lo sites 
is repeated throughout the lattice. In view of the boro- 
carbide structures, we will set La = 1, and let Lq = 1 or 
2 throughout this paper; also, from now on, all energies 
will be measured in units of the hopping t. 

Superconductivity within this model then arises due 
to the balance between formation of pairs and their co- 
herent hopping controlled by e; e.g., as £ — » oo, the sys- 
tem aproaches an insulating state, since hopping will be 
hindered by the high energetic cost of occupying the at- 
tractive sites. Also, in the present context we anticipate 
that the strength of the attraction, \U\, tracks the Debye 
temperature, Gu, in the sense that for a given transition 
metal it increases as the ionic radius of R increases 

The ground state of the Hamiltonian has been 
obtained through Lanczos diagonalization (l7| on lat- 
tices with Ns sites and Nf, electrons (thus the density 
p — Ne/Ns); we take closed chains with a phase change 
(or magnetic flux through the ring), 0, imposed on the 
hopping term between sites Ng and 1, in order to select 
closed-shell configurations i '^sll as to calculate the 
Drude weight (see below). Depending on the SL config- 
uration (Lq =1 or 2), and filling factor, we were able to 
reach lattices as large as Ng = 24; note, however, that 
not all configurations fit into all sizes and occupations 
considered. 

For a given SL configuration, the nature of the ground 
state is determined through the calculation of the follow- 
ing quantities: 

1) The charge gap, which probes single-particle excita- 
tions, defined as 

A = E{Ns,N, + 1) + EiNg, N,-l)- 2E{Ns, N,), (2) 

where E{Ns,N') is the ground state energy for a chain 
with Ns sites and N' electrons. 

2) The Drude weight, which measures the strength of the 
DC conductivity peak, defined as 
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2 d(h^ 



(3) 



= 00 



where, for a given SL configuration, Eq is the ground 
state energy for a system with iVe electrons and TVs sites; 
00 is the phase value which minimizes Eq. 

For both A and D we obtain a series of values for dif- 
ferent system sizes and extrapolate towards Ng — > oo. A 
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Figure 1: Extrapolated Drude weight (left vertical axis) and 
charge gap (right vertical axis) as functions of the site energy 
for the Hubbard superlattice with Lq = La — 1, U = —4, 
and an electronic density p = 5/3. The vertical dashed line 
separates the metallic (M) and superconducting (S) phases. 



vanishing limit of A is associated with a metallic state, 
whereas a non-zero value is associated with either an in- 
sulating state or with a superconducting state; indeed, 
in the case of the homogeneous (i.e., non-layered) attrac- 
tive Hubbard model, it has been exactly shown that the 
single-particle excitation spectrum has a gap for arbi- 
trary densities ■ ^ vanishing limit of D is associated 
with an insulating state, whereas a non-zero value is as- 
sociated with a conducting state, i.e., metallic or super- 
conductor. Therefore, the nature of the ground state is 
determined by the limiting behaviour of both A and D, 
as summarized in Table 

3) The s-wave .singlet (SS) pairing correlation function, 
defined as 



H.C.), 



(4) 



which is used as a consistency check for a superconduct- 
ing ground state. Typically, C{i;£) should display a much 
slower spatial decay in the superconducting state than in 
the normal state. 

Our results are presented by first fixing a value of 
U — —4 (in the intermediate coupling regime), and com- 
paring the data for different values of e and of the elec- 
tronic density; we then let U vary. In Fig. ^ we show 
the extrapolated Drude weight and charge gap for a SL 
with Lo = La ~ 1, and electronic density p — 5/3. Ac- 
cording to Table ^ a phase transition between a metal 
(M) and a superconductor (S) takes place at Ec — 2. As 
expected. Fig. ^ shows that as e is further increased, the 
system approaches an insulating state. Since the insulat- 
ing (I) state is not relevant to the borocarbides (they are 
metallic in their normal state), we will not discuss the 
S-I transition any further. 
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Figure 2: Phase diagram for fixed U — —4, and for a free 
layer with (a) one and (b) two sites. Data points for £c(p), and 
associated error bars, come from gaussian fits to d^A/de^. S, 
M, and I stand for superconducting, metallic, and insulating 
phases, respectively. Full curves guide the eye, and vertical 
dashed lines locate the insulating behavior at pc and e ^ Ec- 
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Figure 3; Log-log plot of the spatial decay of the pairing cor- 
relation function, normalized at ^ = 2, for a 24-site lattice and 
p — 11/6: dashed, dotted, and chain curves represent data for 
homogeneous systems with Uh = —4, 0, and 4, respectively; 
for an SL with Lo = 1 and U = —4, the circles and squares 
are data for e = and 4, respectively. 



Still for U = —4, we calculated the Drude weight and 
the charge gap for other band fillings, both for Lq — I 
3/2, 4/3, 7/4 and 11/6), and for Lq = 2 {p =4/3, 
5/3, and 11/6); the analysis yields the critical boundaries 
edp) depicted in Fig.|21 Several aspects of the phase di- 
agrams are worth stressing. First, a non-zero site energy 
is needed to stabilize the superconducting state. Indeed, 
given the local (i.e., on-site) character of the attraction, 
the presence of one or two free sites in between attrac- 
tive ones would normally (i.e., for s — 0) suppress co- 
herence between pairs. A repulsive site energy on the 
attractive sites then has the effect of spreading the pair 
wave function out to neighboring sites, increasing their 
overlap: pair coherence is thus recovered. Secondly, a 
superconducting ground state is only stable above a crit- 
ical density, pc, the value of which can be obtained from 
a strong couphng (i.e., e :s> \U\ :s> 1) analysis: as one fills 
up an empty lattice, electrons are first accommodated on 
the free sites, so their interaction is only effective when 
the attractive sites start being occupied, i.e., for densi- 
ties above pc = 2Lq/{Lq + La)- Thirdly, at pc the Drude 
weight (charge gap) seems to vanish (to be finite) for all 
£, except at Sc, where it is non-zero (zero); at Ec the 
system is particle-hole symmetric and charge is evenly 
distributed throughout the lattice pd| . 

We now illustrate the analysis of the pairing correla- 
tion function. For a clearer comparison of the different 
rates of decay, in Fig. |3| we eliminate oscillations due 
to the SL structure by plotting the correlations at ev- 
ery other site; also, we normalize the functions such that 
the correlation between an attractive site and its first 
attractive neighbour is set to one; i.e., C(i,£ = 2) = 1, 



with i chosen to be an attractive site. Let us first con- 
sider a homogeneous lattice with — 24 and p = 11/6. 
When Uh = —4 the system is known to be supercon- 
ducting, and, accordingly, the correlation function decays 
about ten times slower than those for the corresponding 
non-superconducting cases, free {Uh = 0), and repulsive 
{Uh =4). Considering now the case of an Lq = 1 super- 
lattice, also with Ng — 24, p = 11/6, and U ~ —4 we see 
that for e = 4, the correlations decay exactly as in the 
corresponding homogeneous system, whereas for e = 
the rate of decay follows more closely that of the free 
system. This confirms the prediction from the analysis 
of A and D, which places the former on the supercon- 
ducting side of the diagram and the latter on the metallic 
side (see Fig. 12)). 

The above procedures are repeated to determine the 
behavior with U for fixed densities. The resulting phase 
diagrams, ec{U), for Lq — 1 are shown in Fig. ^for p = 
5/3 and 3/2. For \U\ > 2, the site energy needed to 
stabilize the superconducting phase increases with \U\, 
due to the spreading effect mentioned earlier. As \U\ — > 
0, a sharp rise of Ec is expected, since one should have a 
metallic state for all e finite when U = Q. For Lq = 2, 
the M-S boundary is qualitatively similar, but the S-I 
boundary is more strongly affected by the presence of a 
second free layer, manifested by a shift to much lower 
values of e. 

A qualitative connection with experiments on borocar- 
bides RT2B2C can now be established, based on the as- 
sociation of |J7| and e with measurable parameters such 
as the Debye temperature and the band offset. Since 
0^1 sets the scale of the BCS zero temperature gap, 
a reasonable starting point is to order the Rs accord- 
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Figure 4: Constant-density sections of the phase diagram 
(e,|?7|,p), for the superlattice with Lo = 1 and (a) p — 5/3, 
and (b) p = 3/2. The full curves guide the eyes through 
the boundary between superconducting (SUC) and metallic 
phases. The chemical symbols correspond to R in _Rr2B2C 
compounds; see the text. 

ing to their Debye temperatures; these, in turn, track 
their ionic radii. Then, in Fig. ^a) we first place the 
whole i?Ni2B2C series in such way that the S-M bound- 
ary lies between Dy and Tb; this pins R to positions on 
the horizontal axis |2lj |. Interestingly, this parametriza- 
tion subsumes the roles of N{Ef), the electron-phonon 
interaction and of O^i, resolving the puzzling fact that 
LaNi2B2C is not superconducting, in spite of the rela- 
tively large N{Ep) and Qd- 

As one changes the transition metal (keeping R fixed), 
by going down a T-column in the periodic table, the elec- 
tronic density in the valence orbitals does not change. 
Therefore, given our model parameters, one expects the 
band offset, e, to be the one most strongly affected by 
this change, since it is related to the energy of the atomic 
levels. Data taken from atomic energy calculations (see, 
for example, indicate that the energy of the higher 
level increases as one goes from Ni to Pd to Pt, as well 
as from Co to Rh to Ir. Accordingly, in Fig. 0fa) the 
horizontal line representing the T = Pt series should be 
placed above the one for Ni, and in such a way that La 
lies within the superconducting phase. This resolves the 
second puzzle of the borocarbides, the appearance of su- 
perconductivity as T is changed: as one moves down a 
column of the VIIIA group, the range of Rs leading to su- 
perconductivity widens. One can therefore predict that 
if a chemically stable phase of, say SmPt2B2C, is grown, 
it should be superconducting. The Pd series suffers from 
a similar scarcity of data due to chemical instabilities: 
apart from superconducting YPd2B2C, only compounds 



with La and Pr have been grown so far, which are not 
superconducting. According to the phase diagram, this 
was already expected, since the horizontal line represent- 
ing the Pd series should lie in between those for Ni and 
Pt. 

We now fix R and move horizontally on the periodic 
table. While the full effects of the change from, say, Ni 
to Co should only be captured within a multi-band treat- 
ment of the transition metals, one expects that the most 
evident feature is a decrease in band filling, possibly ac- 
companied by a renormalization of s. Thus, we place the 
borocarbides with Co, Rh and Ir in the section of the 
phase diagram corresponding to a smaller p, as in Fig. 
2Jb). Since all members of the T = Co series have so far 
failed to display superconductivity, they should be repre- 
sented by a horizontal line lying entirely in the metallic 
phase; see Fig. ^Jb). Many compounds in the Rh series 
(^Rh > £co) have been grown, also without supercon- 
ducting members; this series is therefore represented by 
a line just above that for Co. The phase diagram of 
Fig. ^b) predicts that superconductivity should be more 
likely for the Ir series (eir > ERh), but this series is even 
more severely plagued by chemical instabilities than the 
Pd and Pt ones, and the only compound grown so far 
is the non-superconducting LaIr2B2C. For completeness, 
one should mention that the above reasoning indicates 
that superconducting YRU2B2C should be placed in a di- 
agram corresponding to p < 3/2. Finally, since the dou- 
ble layered materials seem to be non-superconducting, 
they should all lie below the Sc hne in the corresponding 
Lq — 2 diagram. 

In summary, we have proposed a superlattice model to 
describe superconducting layered materials. The model 
is parametrized in terms of band offset e, on-site at- 
traction U, electron density p, and number of free lay- 
ers Lq. The numerical analysis of the one-dimensional 
case (at zero temperature and without local moments) 
established that superconductivity is possible only above 
a critical density Pc{Lo) and for e > edp, U). The model 
captures the essencial physics of superconductivity in the 
borocarbides family, as evidenced by the available exper- 
imental information [regarding the nature of the ground 
state (i.e., metallic or superconducting)] being consis- 
tently accommodated in projections of the phase dia- 
gram. Consequently, one is able to predict which com- 
pounds should become superconductors, once the barri- 
ers of chemical instability are broken. The interplay with 
magnetic ordering is currently being investigated through 
the addition of a Kondo-like term, which couples the con- 
duction electrons to the local moments. 
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